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I Zw 18 is a prototype Blue Compact Dwarf (BCD), characterized by a strong starburst and extremely low metallicity (Z ~ 0.02 Z©). 
It has long been considered a candidate young galaxy in the Local Universe, but recent studies indicate the presence of old stars. 
We analysed archival VLA observations of the 21 cm line and found that the HI associated to the starburst region forms a compact 
fast-rotating disk. The HI column densities are very high, up to ~50-100 M© pc~^ (~ 0.6 - 1.2 X 10^^ atoms cm~^). The rotation 
curve is flat with a steep rise in the inner parts, indicating the presence of a strong central concentration of mass. Mass models with a 
dark matter halo show that baryons may dominate the gravitational potential in the inner regions. A radial inflow/outflow motion of 
~15 km s~^ is also present. 

I Zw 18 appears structurally different from typical dwarf irregulars in terms of gas distribution, stellar distribution and dynamics. It 
may be considered as a "miniature" high-surface-brightness disk galaxy. These dynamical properties must be tightly related to the 
starburst. They also shed new light on the question of the descendants of BCDs. 

There is also extended HI emission towards the outlying stellar complex I Zw 18 C and a ~ 13.5 kpc HI tail. An interaction/merger 
between gas-rich dwarfs is the most likely explanation for the starburst. 
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1. Introduction 

Blue compact dwarfs (BCDs) are low-mass galaxies that are ex- 
periencing a starburst. They are usually characterized by small 
physical sizes (~ 2-3 kpc), low metallicities (0.2 < Z/Z© < 0.02) 
and relatively large amounts of gas (Mhi/Lb > 1). The ques- 
tion was raised whether they are young galaxies undergoing 
their first burst of star formation (Searle & Sargent 1972), but 
several studies based on surface brightness and color profiles 
(e.g. Gil de Paz & Madore 2005) as well as color-magnitude 
diagrams of resolved stellar populations (e.g. Tosi 2009) have 
demonstrated that BCDs contain also old stars, with ages >2- 
3 Gyr. The star-formation histories of the nearby BCDs, as de- 
rived using color-magnitude diagrams (e.g. Tosi 2009), show 
that the starburst is a short-lived phenomenon, typically sus- 
tained for -100 Myr (e.g. McQuinn et al. 2010). Thus, BCDs are 
transition-type dwarfs, but it is not clear whether there are evo- 
lutionary connections with dwarf Irregulars (dirrs), Spheroidals 
(dSphs) and/or Ellipticals (dEs) (e.g. Papaderos et al. 1996; van 
Zee et al. 2001). Also, the mechanisms that trigger, sustain and 
quench the starburst activity are not understood. 

Regarding the H I distribution and kinematics, various stud- 
ies (e.g. van Zee et al. 1998a, 2001) highlighted two striking 
properties of BCDs: i) they show strong concentrations of H I to 
the starburst region near the galaxy centre; ii) they usually have 
steep central velocity gradients. Both properties are not observed 
in more quiescent dIrrs (e.g. S waters et al. 2002). This suggests a 
close connection between the starburst, the compact distribution 
of baryons (gas and stars) and the HI kinematics. The nature of 
the steep velocity gradients is not clear and two main interpre- 
tations are possible: i) fast rotation (e.g. van Zee et al. 2001), 
ii) gaseous inflows/outflows (e.g. Kobulnicky & Skillman 2008, 



Cannon et al. 2004). Fast rotation in the inner regions would in- 
dicate the presence of a strong concentration of mass, that may 
be either luminous or dark. Gaseous inflows/outflows would be 
linked to fuelling processes and/or feedback mechanisms. 

On larger scales, BCDs usually show extended and diff'use 
HI structures, that may form reservoirs for fueling the star- 
burst. Generally, two diff'erent kinds of structures are observed: 
i) extended HI disks in regular rotation, e.g. NGC 2915 (Elson 
et al. 2010), NGC 2366 (Oh et al. 2008); ii) complex filamen- 
tary structures, e.g. II Zw 40 (van Zee et al. 1998a), NGC 5253 
(Kobulnicky & Skillman 2008). The study of these extended HI 
structures can provide key information on the triggering mech- 
anism (external vs internal processes), the properties of the pro- 
genitor galaxies (gas-rich dIrrs vs gas-poor dEs/dSphs), and the 
possible presence of massive gas inflows/outflows. 

We present here a H I study of I Zw 18, the BCD prototype 
(e.g. Zwicky 1966; Searle & Sargent 1972). I Zw 18 is one of 
the most metal-poor galaxies known (12-hlog(0/H)~7.2, Izotov 
& Thuan 1999) and has long been considered a candidate young 
galaxy in the Local Universe, formed within the last 0.5 Gyr 
(e.g. Papaderos et al. 2002; Izotov & Thuan 2004). However, 
Aloisi et al. (2007), using HST observations, detected stars older 
than 1-2 Gyr and ruled out the possibility that I Zw 18 is a 
truly primordial galaxy (as also suggested by e.g. Aloisi et al. 
1999, Ostlin & Mouhcine 2005). The intense star-formation ac- 
tivity started only -20 Myr ago (e.g. Aloisi et al. 1999). Two key 
questions remain open: What triggered the starburst? Why is the 
metallicity so low? 

Previous HI studies (Viallefond et al. 1987; van Zee et al. 
1998b) showed that I Zw 18 is characterized by a strong central 
concentration of H I and a steep velocity gradient, as is typical 



1 



Lelli et aL: Dynamics of Starbursting Dwarf Galaxies: I Zw 18. 
Table 1 : VLA observing parameters 



Project Code Array Conf. Observing Dates Time on Source Calibrators 

AC0710 A 8, 9, 14, 18, 19, 28 Nov & 4, 6 Dec. 2004 333h 0834+555, 1331+305, 0542+498 

AZ0074 B 26 Oct. & 3, 15 Dec. 1995 15.5 h 0834+555, 1331+305, 0542+498, 0137+331 

AP264 C 27 Aug. 1993 6.6 h 0834+555,1331+305 

AP264 D 22 Dec. 1993 2.3 h 0834+555,1331+305 



Table 2: HI Datacubes 



Cube 


Robust 


UV Taper 


Pixel size 


Synthesized beam 


Beam PA. 


Smoothed beam 


AV 


Rms noise 






(U) 


(arcsecxarcsec) 


(arcsecxarcsec) 


(degrees) 


(arcsecxarcsec) 


(km s-i) 


(mJy/beam) 


Low Res. 





10 


3x3 


14.6 X 14.4 


32.7 


20x20 


5.2 


0.40 


Intermediate Res. 





60 


0.5 X 0.5 


3.3x2.9 


89.3 


5x5 


5.2 


0.26 


High Res. 


-1 


150 


0.25 X 0.25 


1.5 X 1.4 


-76.3 


2x2 


10.4 


0.16 



for BCDs. Also, the optical galaxy is surrounded by an extended 
HI envelope that was described by van Zee et al. (1998b) as "a 
fragmenting H I cloud in the early stages of galaxy evolution". 

We analysed archival HI data in order to: i) investigate the 
relation between the gas distribution and the starburst activity; ii) 
clarify the nature of the velocity gradient; iii) study the structure, 
kinematics and origin of the extended gas. 



2. Data reduction & analysis 

We analysed public H I data taken from the VLA archive. The 
observations were carried out between 1993 and 2004, using the 
VLA in all 4 configurations (see table 1). Data from the B, C, 
and D configurations were presented by van Zee et al. (1998b). 
In this new analysis, we included also data taken in 2004 with the 
high-resolution A-array configuration. The correlator was used 
in 2AD mode, with a total bandwidth of 0.8 MHz (-165 km s"^). 
An on-line Hanning taper was applied to the data, producing 127 
spectral line channels with a width of 6.3 kHz (~1.3 km s"^). 

The raw UV data were interactively flagged, calibrated and 
combined using the AIPS package and following standard VLA 
procedures. Next, the UV data were mapped using a robust 
weighting technique (Briggs 1995) and various Gaussian base- 
line tapers to attenuate the longest baselines. We built three da- 
tacubes with difl'erent spatial resolutions and pixel sizes by us- 
ing difl'erent combinations of the Robust parameter and the taper 
FWHM (see table 2). After various trials, we chose the combi- 
nations that minimize sidelobes and wings in the beam profiles. 

After the Fourier transform, the datacubes were analysed 
using the Groningen Imaging Processing SYstem (GIPSY) 
(van der Hulst et al. 1992). Continuum maps were constructed 
by averaging line-free channels. Because of the small bandwidth 
of the observations, few line-free channels were available and 
the resulting continuum- subtracted datacubes showed correlated 
noise in the spectral direction. Thus, we constructed a contin- 
uum map by using a mask, defining the area of H I emission in 
every channel and averaging, for each spatial pixel, all the chan- 
nels without H I signal. The masks were constructed by smooth- 
ing the datacubes both in velocity (by a factor 4) and spatially 
(by a factor ~3, i.e. at 45'', 10'' and 5" for the low, intermedi- 
ate and high resolution data, respectively) and clipping at 2.5(T^ 
(where cTs is the noise in the smoothed cubes). The masks were 
inspected channel by channel and remaining noise peaks were 
blotted out. 

The use of a mask for the continuum subtraction may have 
the disadvantage that the noise is no longer uniform across the 
channel maps, as a diff'erent number of channels is used at every 



pixel to build the continuum map. Thus, we built signal-to-noise 
(S/N) maps for every channel (similarly to Verheijen & Sancisi 
2001) and calculated a pseudo-1 cr contour by averaging the va- 
lues of the pixels with 0.75<S/N<1.25. The resulting pseudo-lcr 
level is close to that obtained by calculating the noise in a box 
without signal, suggesting that the noise is still almost uniform. 

The channel maps were cleaned (Hogbom 1974) down to 
0.5(T, using the masks to define the search areas for the clean- 
components, which were then restored with a Gaussian beam of 
the same FWHM as the antenna pattern. Next, to boost the S/N- 
ratio, the cubes were smoothed in velocity to a resolution of 5.2 
km s"^ (10.4 km s"^ for the high-resolution data) and spatially to 
20", 5" and 2" for the low, intermediate and high resolution data, 
respectively. Table 2 summarizes the properties of the cubes. 

Total H I maps were constructed by summing the signal in- 
side the clean-masks. A pseudo-3(T contour was calculated fol- 
lowing Verheijen & Sancisi (2001). Velocity fields were built 
by fitting a Gaussian function to the HI line profiles. Fitted 
Gaussians with a peak intensity less than 2.5cr and a FWHM 
smaller than 5.2 km s"^ were discarded; remaining noise in 
the velocity fields (i.e. signal outside the pseudo-3(T contour of 
the total HI maps) was blotted out. The HI line profiles are 
quite broad and asymmetric, thus the velocity fields must be 
considered just as a rough indication of the global kinematics. 
Our kinematical analysis is based on Position- Velocity diagrams 
(Sect. 3) and on 3D models of the observations (Sect. 4.1). 



3. Hi distribution and (cinematics 

In this section we describe the overall HI structure of I Zw 18. 
We adopt the standard nomenclature introduced by Davidson 
et al. (1989). The main body is designated as I Zw 18 A (Fig. 
2, top-right) and is characterized by two starburst regions: one to 
the North- West (NW) and one to the South-East (SE). The fight 
concentrations denoted by Davidson et al. (1989) with B, D and 
E are background galaxies. The stellar complex to the NW is 
named I Zw 18 C or C-component. We assume a distance of 18.2 
Mpc, as derived from the tip of the red giant branch (Aloisi et al. 
2007) and confirmed by observations of Cepheids (Fiorentino 
et al. 2010). 

We use data at three diff'erent resolutions (see table 2) in or- 
der to probe diff'erent spatial scales and HI column densities. 
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Fig. 1: Channel maps at a resolution of 20'' (top) and 5'' (bottom). Red-solid contours are at 1.5, 3, 6, 12, 24 x cr. Red-dashed 
contours are at -3, -1.5 x cr. Black contours show two isophotes of a B-band image (from Gil de Paz et al. 2003); the object to the 
North-East (R.A. = 9h 30m 34s, DEC = 55° 28'') is a foreground star. 
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Fig. 2: Top: integrated H I maps at a resolution of 20'' (left) and 5'' (right), overlayed on a HST image (from Aloisi et al. 2007). The 
box in the left panel shows the area covered by the right panel. The dashed line corresponds to the pseudo-1. Scr density contour. In 
the map at 20'', contours are at 0.25 (dashed), 0.5, 1,2, 4, 8, 16x 10^^ atoms cm"^. In the map at 5", contours are at 3 (dashed), 
6, 12, 24, 48 X 10^^ atoms cm"^. The circle shows the beam size. Middle: velocity fields at a resolution of 20" (left) and 5" 
(right). The box in the left panel shows the area covered by the right panel. Contours range from 722.4 to 805.6 km s"\ with steps 
of 5.2 km s"^ The circle shows the beam size. The dashed line shows the path followed to obtain the position- velocity diagram. 
Bottom: position- velocity diagrams at a resolution of 20" (left) and 5" (right). Contours are at -1.5 (dashed), 1.5, 3, 6, 12, 24 x cr. 
Jhe box in the left panel shows the region covered by the right panel. The vertical line corresponds to the cross in the velocity fields. 
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R.A. (B1950) R.A. (B1950) R.A. (B1950) 

Fig. 3: Left: integrated HI map at 2" resolution. Countours are at 3 (dashed), 6, 9, 12, 15 x 10^^ atoms cm"^. Middle: Ha image 
(from Cannon et al. 2002) overlayed with the integrated H I map at 1" resolution. Right: velocity field at 1" resolution. Contours 
range from 722.4 to 805.6 km s"^ with steps of 10.4 km s"^ The circle shows the beam size. 



3. t. The low-resolution view 

The low-resolution data (with FWHM = 20'' -1.8 kpc and 3cr 
column density sensitivity Nhi(3(t) ~ 4 x 10^^ atoms cm"^ per 
channel) illustrate the large-scale overall structure of I Zw 18. 

Figure 1 (top) shows the channel maps, overlayed with two 
isophotes of a B-band image (from Gil de Paz et al. 2003). The 
central HI emission presents a velocity gradient at a position 
angle P.A.^140°-150°. At velocities from -750 to -720 km s"^ 
extended emission appears also to the South. 

Figure 2 (top-left) shows the integrated H I map, overlayed 
on a V-band HST image (from Aloisi et al. 2007). I Zw 18 A is 
associated with a strong concentration of gas, while diff'use emis- 
sion extends beyond the optical galaxy, covering an angular size 
of ~3'.5 (-18.5 kpc). The HI gas to the South of I Zw 18 A dis- 
plays a tail-like morphology that extends over ~ 2'. 5 (-13. 5 kpc). 

Figure 2 (middle-left) shows the velocity field. The main 
body (I Zw 1 8 A) is associated with the central velocity gradient. 
The southern "tail" does not seem to be kinematically connected 
with the SE region of I Zw 1 8 A, as the gas velocity changes 
abruptly from -790 km s"^ to -720 km s"^ Moreover, at the 
junction between I Zw 18 A and the "tail", the HI line profiles 
are double peaked, suggesting that there are two distinct compo- 
nents, possibly well- separated in space but projected to the same 
location on the sky. 

Figure 2 (bottom-left) shows a Position- Velocity (PV) dia- 
gram, obtained from the 20'' datacube following the "tail" (the 
dashed line overlayed on the velocity field). The central velocity 
gradient (associated with I Zw 1 8 A) is very steep and there is 
a spatial broadening towards the NW direction between -700 
and -780 km s"^ The gas to the South forms a coherent kine- 
matical structure at velocities between -710 and -760 km s"^ 
Strikingly, the broadened part of the PV-diagram and the "tail" 
are almost at the same velocities, suggesting that they may be 
physically connected. This hypothesis will be investigated fur- 
ther in Sect. 5. 

3.2. The intermediate-resolution view 

The intermediate-resolution data (FWHM = 5" - 440 pc and 
Nhi(3(t) - 4 X 10^^ atoms cm"^ per channel) illustrate the HI 



emission associated with I Zw 18 A and I Zw 18 C, and their 
possible connections. 

Figure 1 (bottom) shows the channel maps, overlayed with 
two isophotes of a B-band image. The H I emission to the NW is 
spatially resolved but still visible, whereas the southern "tail" is 
completely resolved out, indicating the diff'use nature of this gas. 
Between -770 and -740 km s"^ there are HI clumps near the 
C-component. Their association with I Zw 18 C is very likely, 
because the H I clumps are at the same velocities as the Ha emis- 
sion (ysys,Ha[C] = 751 ± 5 km s-\ Dufour et al. 1996). 

Figure 2 (top-right) shows the integrated H I map, overlayed 
on a V-band HST image. The main body is characterized by two 
H I peaks, roughly corresponding to the NW and SE starburst 
regions. With respect to I Zw 18 A, the HI is more extended to 
the North- West, in the direction of I Zw 18 C. 

Figure 2 (middle-right) shows the velocity field, while fi- 
gure 2 (bottom-right) shows a PV-diagram obtained following 
the dashed line on the velocity field. The steep velocity gradi- 
ent is approximately along the two HI peaks. The HI emission 
to the NW shows a shallow velocity gradient from I Zw 18 A 
to I Zw 18 C and seems to connect the two stellar bodies. The 
presence of a "gaseous bridge" is supported also by Ha obser- 
vations: Dufour & Hester (1990) detected diff'use Ha emission 
connecting I Zw 18 A and I Zw 18 C, while Dufour et al. (1996) 
measured Ha velocities between the two bodies that are very 
similar to those observed in H I (their Fig. 3). Along I Zw 18 A, 
instead, the Ha velocity gradient shows a "wiggly" behaviour 
that is not observed in H I. This may be due to the presence of an 
Ha superbubble (see Martin 1996 and Sect. 5.2). Also, there is 
H I emission to the West of the main body, that shows a velocity 
gradient and seems to have an Ha counterpart (see Sect. 5.2). 

3.3. The high-resolution view 

The high-resolution data (FWHM = 2" - 180pc andNni (3o-) - 
2 X 10^^ atoms cm"^ per channel) show in detail the H I emission 
associated to the NW and SE starburst regions. 

Figure 3 shows the total H I map in grayscale (left) and in 
contours overlayed on a Ha image (middle). The two HI peaks 
are spatially resolved. The H I clump to the SE coincides with 
a complex of H II regions and has a strong peak, where the H I 
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RA (B1950) RA (B1950) 

Fig. 4: Left: velocity field at 5'' resolution. The ellipse shows the modelled disk. The centre, the major and the minor axis are shown. 
The box shows the area covered by the right panel. The circle shows the beam. Countours are the same as in Fig. 2. Right: HST 
image overlay ed with the observed H I map at 2" (black) and the H I map obtained from a 3D model with / = 70° and zo = 100 pc 
(red). Countours are at 3 (dashed) and 6 (solid) x 10^^ atoms cm"^. The cross marks the centre. The circle shows the beam. 



column density reaches -100 M© pc"^ (~ 1.2 x 10^^ atoms cm^). 
In between the two clumps, an HI hole is associated with a 
strong H II region, suggesting that the neutral gas has been con- 
sumed, ionized and/or blown out by young stars. The H I clump 
to the NW coincides with a Hof shell, that surrounds the bulk of 
the young stars (see Fig. 1 of Cannon et al. 2002). This shell is 
probably connected with the high- velocity Ha emission detected 
by Dufour et al. (1996) and Martin (1996) at ±200 km s'^ with 
respect to the systemic velocity. 

The velocity field at 2" resolution is shown in Fig. 3 (right). 
This velocity field is very uncertain because of the clumpy H I 
distribution, the asymmetric line profiles and the low S/N-ratio 
of the data at this high angular resolution. However, it shows a 
clear velocity gradient from the SE to the NW region, as already 
observed at lower resolutions. The H I depression is in the ap- 
proaching NW side of the galaxy. 

4. Dynamics of I Zw 18 A 

In Sect. 3 we described the overall structure of I Zw 18. Two 
main facts need to be explained: i) the nature of the steep velocity 
gradient associated with I Zw 18 A; ii) the origin of the extended 
H I emission to the South and to the North- West of the galaxy. In 
this section we will focus on the dynamics of I Zw 1 8 A, while 
in Sect. 5 we will study the large-scale gas emission. 

4A. Kinematical models 

There is some controversy in the literature regarding the velocity 
gradient of I Zw 18 A: Viallefond et al. (1987), Petrosian et al. 
(1997) and van Zee et al. (1998b) analysed 2D velocity fields 
and interpreted the gradient as rotation, whereas Skillman & 
Kennicutt (1993) and Dufour et al. (1996) obtained ID long-slit 
spectroscopy and argued that the gradient may result from the 
merger of two (or more) gaseous clouds. The velocity gradient 
is along the optical major axis of the galaxy (Fig. 4) and the velo- 



city field shows the typical pattern due to rotation: this strongly 
suggests the presence of a rotating disk. Here we present 3D 
kinematical models, which demonstrate that the H I disk is dif- 
ferentially rotating and has a global inflow/outflow motion. 

The disk is modelled by a set of gas rings, with fixed cen- 
tre, systemic velocity, position angle, inclination, surface den- 
sity, thickness, velocity dispersion and rotation velocity. The 
centre, the systemic velocity and the position angle were esti- 
mated by eye using both optical and HI data (see table 3). The 
centre is between the NW and the SE starburst regions (see Fig. 
4, right). For the radial distribution, we used the H I surface den- 
sity profile, derived from the total H I map at resolution by 
azimuthally averaging over ellipses (Fig. 6, top). For the verti- 
cal distribution, we assumed an exponential law exp(-z/zo)- We 
built a set of models with diff'erent values for the inclination /, the 
scale height zo and the velocity dispersion (Thi, assuming that 
each of these parameters is constant with radius. The inclination 
and the scale height are constrained by the observed HI map; 
their values are slightly degenerated but do not strongly aff'ect 
the final result; we assumed / = 70° and zo = 100 pc (see Fig. 4, 
right). The mean velocity dispersion is constrained by the shape 
of diff'erent PV-diagrams and values higher than ~10 km s"^ are 
ruled out; we assumed crm =7.5 km s"^ . A thickness of 100 pc 
and a H I velocity dispersion of 7.5 km s"^ are typical values for 
a HI disk. 

The actual H I distribution of I Zw 1 8 A is clearly not axi- 
symmetric (see Fig. 3). Thus, once the structural and geomet- 
rical parameters of the disk were fixed, we built models with 
axi-symmetric kinematics but a clumpy H I distribution, i.e. the 
surface density varies with position as in the observed H I map. 
The procedure is as follows: we built a disk model with uniform 
density distribution (fixing zo, crm and the rotation curve), pro- 
jected it on the sky and then renormalized the H I line profiles at 
every spatial position in order to reproduce the H I density distri- 
bution observed at resolution. For the rotation curve, we tried 
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Fig. 5: Comparison between different 3D kinematical models and the observations. Top: Position- Velocity diagrams at a resolution 
of 2" and 5''. The slice position are the major and minor axis, as indicated in Fig. 4 (left). Contours are at -1.5 (dashed), 1.5, 3, 
6, 12 X (T. Bottom: Channels maps at 1" resolution (red contours) overlayed with two models: differential rotation (black-solid 
contours) and solid body rotation (black-dashed contours). Contours are at 2(T. The cross marks the galaxy centre. The circle shows 
the beam size. See text for details. 



7 



Lelli et aL: Dynamics of Starbursting Dwarf Galaxies: I Zw 18. 



Table 3: Properties of I Zw 18. 



Parameter 


T 'Vwi ISA 


T 'Vwi 1 8 




JU jU. J + U. i 


rvnhorvnio'7 so rv o 
JU 11 . ^ ± U.Z 




jj LI ^1 + i 


JJ zo Uo + z 


»^sys V^Al 5> ) 


757 + 4 


751+5 


Position Angle (°) 


145 + 5 




Inclination Angle (°) 


70 + 4 




Vrot (km s-i) 


38 + 4.4 




U (10^ Lo) 


13.9 


1.1 


Lr (10^ Lo) 


5.8 


0.4 


Mhi (10^ Mo) 


1.0 


<1.2 


Mdyn (10^ Mo) 


3 + 1 





Notes. Optical luminosities were calculated using the apparent magni- 
tudes from Papaderos et al. (2002), the distance from Aloisi et al. (2007) 
and the solar absolute magnitudes from Binney & Merrifield (1998). 
The HI mass of I Zw 18 C refers to the extended emission described 
in Sect. 5, althought only part of this gas may be physically associated 
with the C-component. The southern "tail" accounts for ~ 0.5 x 10^ Mo. 

two extreme cases: solid body (slowly rising rotation curve) and 
differential (steeply rising and flat rotation curve). 

Figure 5 (top) shows PV-diagrams obtained from difl'erent 
models and from the observations at a resolution of and 5''. 
The slice positions are the major and minor axis, as indicated in 
Fig. 4 (left). The velocity gradient along the major axis is grossly 
reproduced by all the models. This demonstrates that: i) a rotat- 
ing disk is a good representation of the data; ii) the asymme- 
try between the NW and the SE region is mostly caused by the 
clumpy H I distribution. Moreover, it is possible to discriminate 
between solid body and difl'erential rotation. The observed PVD 
along the major axis shows HI emission close to the galaxy 
centre {R < 5'') at high rotational velocities (-810 and -730 
km s"^). The solid body model does not reproduce such emis- 
sion, as the gas in the inner radii is mostly concentrated near the 
systemic velocity. The difl'erentially rotating model, instead, cor- 
rectly reproduces the high- velocity gas. This is clearly illustrated 
by the channel maps at resolution at receding velocities (Fig. 
5, bottom): the solid body model (dashed line) is not extended 
enough towards the galaxy centre, whereas the difl'erentially ro- 
tating model (solid line) gives a good match with the observa- 
tions. The approaching NW side of the galaxy is not reproduced 
as well as the receding SE side. The NW starburst region is more 
active than the SE one (Contreras Ramos et al. 2011) and this 
may afl'ect its local kinematics. Nonetheless, also on this side, a 
difl'erentially rotating disk is preferable to a solid body one. 

A simple rotating disk (independently from the assumed ro- 
tation curve) cannot reproduce the observed PV-diagrams along 
the minor axis (Fig. 5, top), because of the kinematic asymmetry 
and the presence of H I emission at velocities forbidden by cir- 
cular motions. This efl'ect is visible also in the velocity field at 5'' 
resolution (Fig. 4, left): the kinematic minor axis, defined by the 
contours close to the systemic velocity, is not orthogonal to the 
kinematic major axis. Usually, this is attributed to radial motions 
(e.g. Fraternali et al. 2002). Alternatively, the non-orthogonality 
between the minor and major axes may indicate the presence of 
a bar-like or oval distortion (e.g. Bosma 1978). The optical im- 
ages of I Zw 18, however, do not indicate the presence of such 
strong distortions. We improved the diff'erentially rotating model 
by adding a global radial motion of 15 km s"^. The resulting 
model reproduces the HI emission at forbidden velocities. The 
model, however, cannot reproduce all the details present in the 
observed PV-diagram. The non-circular motions are not uniform 
across the disk and small variations (of the order of 3-4 km s"^) 



could account for the observed discrepancies. It is not possible 
to discriminate between inflow and outflow, as it is not known 
which side of the disk is the near one. Vertical motions with 
roughly the same speed as the radial ones may also be present. 
Non-circular motions in excess of 20 km s"^ are ruled out, con- 
firming that the disk kinematics is dominated by rotation. 

The rotation curve used to build the best model was not 
derived from a standard tilted-ring fit of the velocity field 
(Begeman 1987), but it was derived with a trial and error ap- 
proach building a series of 3D models. The uncertainties on the 
rotation velocities are difficult to quantify. We made a conser- 
vative estimate of the uncertainties assigning an error equal to 
AV/2.35 = 4.4 km s"^ where AV is the velocity resolution of 
the datacube. The first point of the rotation curve is the most 
uncertain, as the velocity dispersion in the inner ring may be 
higher than the mean value of 7.5 km s"^ For example, if (Thi 
is assumed to be 10 km s"^ higher than the mean value, the ro- 
tation velocity would decrease of -5 km s"^ We calculated the 
asymmetric drift correction following Meurer et al. (1996), but 
it turns out to be smaller than the errors. 

4.2. Mass models 

In Sect. 4.1 we showed that I Zw 18 A has a rotating HI disk. 
The rotation curve is uncertain, but it shows an inner steep rise 
and an outer flat part. This indicates the presence of a strong cen- 
tral concentration of mass, that may be luminous or dark. Using 
this rotation curve, we built mass models to estimate the relative 
contributions of luminous and dark matter to the gravitational 
potential. We followed Begeman (1987). 

The contribution of the gaseous disk was computed using the 
surface density profile derived from the total H I map at reso- 
lution (Fig. 6, top), multiplied by a factor 1 .4 to take into account 
the presence of Helium. The possible gravitational effect of the 
H I components outside the disk (i.e. the "tail" to the South and 
the extensions to the North- West and to the West) was not taken 
into account. Molecular gas was not explicitly considered in the 
mass model because its amount is very uncertain (Leroy et al. 
2007). However, if molecules are distributed as the stars, their 
contribution is reffected in an increase of the stellar mass-to-light 
ratio (M.,/L). Consistently with the models in Sect. 4.1, we as- 
sumed an exponential vertical distribution with zo = 100 pc. 

The contribution of the stars was computed using the R-band 
surface brightness profile from Papaderos et al. (2002) (Fig. 6, 
middle), that is derived from an HST image after the subtrac- 
tion of the nebular emission (dominated by the Ha line). The 
color profiles from Papaderos et al. (2002) (their Fig. 11) show 
that, after the subtraction of the nebular emission, the color of 
I Zw 18 A is almost constant with radius. Thus, it makes sense 
to use a constant value of M^^/Lr. We assumed a stellar disk with 
a vertical density distribution given by p(z) = sech^(z/zo) (van 
der Kruit & Searle 1981), with zo = 100 pc. 

For the dark matter distribution, we assumed a pseudo- 
isothermal halo described by equation: 



where po is the central density and Tc is the core radius, po and 
Tc are free parameters of the mass models. 

Figure 6 (bottom) shows the "maximum disk" decomposi- 
tion of the rotation curve. The resulting stellar mass-to-light ra- 
tio is M=,/Lr ~ 1.5. In the maximum disk hypothesis, the baryons 
dominate the gravitational potential in the inner regions, while 
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Fig. 6: Top: HI surface density profile (inclination corrected), 
derived by ellipse averaging over the entire disk (dots) and over 
the approaching and receding sides separately (down-triangles 
and up-triangles, respectively). Middle: R-band surface bright- 
ness profile (from Papaderos et al. 2002). Bottom: "maximum 
disk" decomposition of the rotation curve. Dots and circles show 
the observed curve. Lines show the contributions due to gas, 
stars and dark matter and the resulting rotation curve. 



the dark matter halo dominates in the outer parts. The param- 
eters of the halo are uncertain because the sampling of the ro- 
tation curve is poor. The halo shown in Fig. 6 (bottom) has: 
po = 833 X 10"^ Mo pc"^ and = 0.13 kpc. 

A MJLr -1.5 implies a stellar mass of -9x10^ M©. 
According to Aloisi et al. (1999), the starburst started ~20 Myr 
ago with a star-formation rate (SFR) of 6x10"^ M© yr"\ giv- 
ing a mass in young stars of -10^ Mq^ Thus, the newly formed 
stars and the concentration of H I cannot explain the steep rise of 
the rotation curve, implying that the mass concentration is due 
to either old stars, or molecules, or dark matter. Old stars have 
been detected (Aloisi et al. 2007) and their total mass can be 
constrained by deriving the galaxy Star-Formation History from 
Color-Magnitude diagrams. The maximum-disk value requires a 
mean SFR of -7x10"^ M© yr"^ over the last 13 Gyr, that can- 
not be ruled out. Regarding molecules, the upper limit for the H2 
mass within -400 pc is ~ 7 x 10^ M© (Leroy et al. 2007)^, using 
a Galactic CO-to-H2 conversion factor (Xco)- However, Leroy 



^ Aloisi et al. (1999) assumed a distance of 10 Mpc. Considering the 
new value of 18.2 Mpc, the starburst parameters slightly change, but the 
mass in young stars remains almost the same (F. Annibali, priv. comm.). 

^ We reseated the original values to a distance of 18.2 Mpc. 
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Fig. 7: Position- Velocity diagrams at a resolution of 20'' (top) 
and 5'' (bottom). They are derived following the path shown in 
Fig. 2. Countours are at -1.5 (dashed), 1.5, 3, 6, 12, 24 cr. The 
star shows the Ha velocity and the spatial position of I Zw 18 C. 



et al. (2007) argued that, in I Zw 18, Xco may be 10"^ times the 
Galactic value. The same result is found by extrapolating the re- 
lation between Xco and metallicity by Boselli et al. (2002) down 
to the metallicity of I Zw 18. Thus, the H2 mass within ~400 pc 
may be dynamically important and as high as -7x10^ M©. 

We also considered Modified Newtonian Dynamics 
(MOND) (Milgrom 1983, Sanders & McGaugh 2002) to 
describe the rotation curve. We fixed ao = 1.21 x 10"^ cm s"^ 
(Begeman et al. 1991) and the distance D = 18.2 Mpc (Aloisi 
et al. 2007), thus the only free parameter is M=,/Lr. MOND 
provides acceptable fits and gives M*/Lr = 1.5 using the "stan- 
dard" interpolation function (Milgrom 1983) and M.,/Lr = 1 
using the "simple" one (Famaey & Binney 2005). 

Following McGaugh (2011), we check the position of 
I Zw 18 A on the baryonic Tully-Fisher relation. The galaxy 
follows the correlation within the observed scatter. 



5. The extended emission 

In this section we study the extended HI emission. This may 
provide some clues to the mechanism that triggered the starburst. 
Also, we compare the large-scale H I and Ha emission to inves- 
tigate the possible presence of outflows. 
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R.A. (B1950) R.A. (B1950) Distance along the slice [arcnnin] 

Fig. 8: The HI emission at 30'' resolution, after the subtraction of the main body. Left: B-band image overlay ed with the total HI 
map. Countours are at 2.4, 4.8, 9.6, 19.2 x 10^^ atoms cm"^. Middle: Velocity field. Contours range from 735.8 to 767 km s"^ with 
steps of 10.4 km s"^ The cross shows the position of I Zw 18 A. The dashed line shows the path followed to obtain the position- 
velocity diagram. Right: Position- Velocity diagram. Countours are at -1.5 (dashed), 1.5, 3, 6, 12 cr, where cr = 0.3 mJy/beam. The 
vertical line corresponds to the cross in the velocity field. The star shows the Ha velocity and the spatial position of I Zw 18 C. 



5.1. The C-component and the Hi tail 

In Sect. 3.1 we reported two puzzling results (see Fig. 2, left): 
i) an HI "tail" at velocities between 710-760 km s"^ extending 
to the South of I Zw 1 8 A and kinematically disconnected from 
the South-East side of the central rotating disk; ii) a broadening 
of the HI emission in the PV-diagram almost at the same ve- 
locities (700-780 km s"^) to the North- West (in the direction of 
I Zw 18 C). This is clearly shown in Fig. 7, where PV-diagrams 
at 20'' and 5" resolution (with diff'erent column density sensitiv- 
ities) are plotted on the same scale. 

In order to study these components and their possible con- 
nection in more detail, we subtracted the compact H I disk of 
I Zw 1 8 A from the surrounding extended H I emission. We used 
the high-resolution data to define the emission from the disk and 
we subtracted this emission from the low-resolution datacube^. 
Subsequently, we smoothed the residual datacube to 30" and to 
10.4 km s"^ and used it to obtain: i) a total H I map by summing 
the channels in the velocity range -700-780 km s"^; ii) a velo- 
city field by estimating an intensity- weighted mean velocity; iii) 
a PV-diagram by following the tail (dashed line in Fig. 8). 

Figure 8 shows the results of the subtraction. Interestingly, 
extended H I emission is centered on I Zw 1 8 C (left panel) and 
forms a coherent kinematical structure with velocities between 
-700 and 800 km s"^ (right panel). The physical association of 
the C-component with this surrounding HI emission is likely 
because the Ha systemic velocity of I Zw 18 C is -751+5 km s"^ 
(Dufour et al. 1996), (see star in Fig. 8, right). 

The southern H I tail seems to be connected in space and in 
velocity with the H I structure around I Zw 18 C, as is shown by 
the velocity field (Fig. 8, middle). The connection may be behind 
or in front of I Zw 1 8 A. Possible interpretations of the extended 
HI emission will be discussed in Sect. 6.1. 



^ Technically, we built a mask containing only the H I signal from 
the disk and we cleaned the datacube at 1.5'' x 1.4'' resolution down 
to Icr, using the mask to define the search areas. We restored the clean- 
components on a blank cube, using a Gaussian beam of 20". The result- 
ing "clean-components" cube contains only the emission from the disk, 
but at the desired resolution of 20". Finally, the "clean-components" 
cube was subtractred channel by channel from the 20" datacube. 
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Fig. 9: Ha image (Gil de Paz et al. 2003) overlayed with the HI 
emission. White contours show the HI emission at 2. ''8 x 3. ''3 
resolution, integrated in the velocity range 800-770 km s"^ and 
correspond to 5, 10, 20, 40 x 10^^ atoms cm"^. The black line 
shows the pseudo-3cr contour of the total H I map at 5'' res. 



5.2. The connection between Hi and Ha emission 

In Sect. 3.3 we pointed out the relative H I and Ha distributions 
in the inner regions of I Zw 1 8 A. Here we compare the distribu- 
tion and kinematics of neutral and ionized gas on larger scales. 

Figure 9 shows an Ha image (from Gil de Paz et al. 2003) 
overlayed with the H I emission (black and white contours). The 
Ha emission extends well beyond the stellar body and is almost 
perpendicular to the H I disk, suggesting the presence of an out- 
flow. This interpretation is in agreement with the Ha long- slit 
spectroscopy by Martin (1996), who detected a bipolar super- 
bubble expanding with velocities of ±60 km s"^ out to ~2 kpc 
from I Zw 1 8 A. The H I to the North-East of the main body (see 
black contour) seems to border the Ha emission and is almost at 
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the same velocities, suggesting that part of the diffuse H I may 
be associated with the outflow. 

The Ha emission presents also a prominent arc to the West 
of I Zw 18 A. This feature was first identified by Dufour & 
Hester (1990) and interpreted as a radiation-bound ionization 
front driven into the ISM. Petrosian et al. (1997), instead, argued 
that the Ha arc contains also stellar emission and suggested that 
it is a structure with stars able to ionize the gas "in situ". The H 1 
to the West of I Zw 18 A is associated with the Ha arc. To show 
this, we summed the channel maps at 2. ''8 x 3/^3 resolution in 
the velocity range 800-770 km s"^ (see Fig. 9, white contours). 
Moreover, the Ha velocity field of Petrosian et al. (1997) shows 
a gradient along the arc similar to the one observed in Fig. 2 
(middle-right), confirming the physical association between Hi 
and Ha. The hypothesis of a radiation-bound ionization front is 
difficult to reconcile with the presence along the Ha arc of high- 
density (-10^^ atoms cm"^) neutral gas. 

6. Discussion 

6.1. Observational evidence & interpretation 
This H I study of I Zw 1 8 has shown that: 

- I Zw 1 8 A has a compact rotating disk with very high H I 
densities. The rotation curve is flat with a steep inner rise, 
indicating the presence of a strong concentration of mass. 
Also, a global inflow/outflow motion is present; 

- I Zw 18 C is located in the direction of the major axis of 
I Zw 18 A and is almost at the same velocities as its ap- 
proaching side. Gas emission with a smooth velocity gradi- 
ent connects the two stellar bodies. I Zw 18 C appears to be 
at the centre of a difl'use H I structure; 

- a H 1 tail extends to the South of I Zw 1 8 A out to ~ 1 3 .5 kpc. 
The tail has a coherent kinematical structure and seems to be 
connected with the H I emission to the North- West. 

Studies of the resolved stellar populations (Aloisi et al. 2007; 
Contreras Ramos et al. 2011) have shown that: i) the two star- 
burst regions in I Zw 18 A (NW and SE) are embedded in a 
common envelope of old stars with ages > 1 Gyr; ii) I Zw 1 8 A 
and I Zw 18 C are two completely separate stellar bodies and 
there are no stars between them; iii) also I Zw 18 C contains 
both old (>1 Gyr) and young (-10 Myr) stars, but its current 
star formation rate (SFR) is lower than that of I Zw 1 8 A. 

For the interpretation, we consider first the hypothesis of 
an interaction/merger of two (or more) gas-rich dwarfs. It is 
well-known that interactions/mergers can produce tidal tails (e.g. 
Toomre & Toomre 1972). Also, numerical simulations (e.g. 
Hibbard & Mihos 1995) suggest that mergers can lead to gas in- 
flows, produce strong gas concentrations and trigger intense star- 
formation. Thus, an interaction/merger may provide an explana- 
tion for: i) the concentration of H I, ii) the on-going starburst, 
and iii) the southern Hi tail. I Zw 18 C may be either a "relic" 
of the interaction or a dwarf galaxy that is interacting/merging 
with I Zw 18 A. The two objects are at a projected distance of 
~2.2 kpc, the difl'erence between their systemic velocities is ~12 
km s"\ and they are connected by HI emission with a smooth 
velocity gradient. The ratio between the R-band luminosities of 
I Zw 18 A and I Zw 18 C is ~ 14, thus this would be classified 
as a minor merger. The merger hypothesis may also explain the 
extremely low metallicity, as discussed in the following. 

Bekki (2008) argued that BCDs with low nebular metallicity 
are the results of mergers between gas-rich dwarfs with extended 
HI disks. According to his simulations, the central starburst is 



fuelled with metal-poor gas transferred from the outer regions 
of the extended disks, where the star formation and the chem- 
ical enrichment were not efficient due to the low HI densities. 
I Zw 18 is consistent with this picture, as the C-component is 
surrounded by an extended H I structure, that does not have a 
stellar counterpart and, likely, has not been efficiently enriched 
by SN explosions. Therefore, this HI structure may provide 
"fresh" unprocessed gas into the starburst regions of I Zw 18 A. 
A similar mechanism of metal dilution has been proposed also 
by Ekta & Chengalur (2010) to explain why I Zw 18 and the 
other extremely metal-deficient BCDs are outliers of the mass- 
metallicity relation. 

The presence of the extended HI emission, including the 
tail, has been considered above as supporting evidence for the 
merger, but could it be instead the result of a blowout from the 
starburst? Indeed, Ha observations of I Zw 1 8 suggest the pres- 
ence of such an outflow (Sect. 5.2). Also, numerical simulations 
predict that starbursting dwarfs undergo massive ouflows be- 
cause they have a shallow gravitational potential (e.g. Mac Low 
& Ferrara 1999). The rate of the outflowing gas dMout/dt can be 
roughly estimated as: 

JMout 2 X £ X SNR X £'sN 



where SNR is the rate of supernovae (SN), £'sn is the mean en- 
ergy of a SN, Vesc is the escape velocity and s is the efficiency 
of the SN-feedback. Thus, Mout = dMp^Jdt x Ar, where M is 
the duration of the starburst. Assuming £^sn = 1.2 x 10^^ erg, 
Vesc = V2 X Vrot, £ = 0.15, AT = 20 Myr (Aloisi et al. 1999) 
and SNR = 0.01 xSFR with SFR = 0.06 M© yr'^ (Aloisi et al. 
1999), we get Mout ~ 5 - 8 x 10^ M©. The total mass of the ex- 
tended gas is ~ 1.6 X 10^ M© (corrected for the presence of He). 
Thus, a massive outflow may explain all the difl'use gas if only 
slightly higher values of S FR, AT and s are assumed. Since the 
extended H I emission is entirely at approaching velocities, any 
outflow should be highly asymmetric and confined. 

Finally, there is also the hypothesis of a "fragmenting HI 
cloud in the early stages of galaxy evolution", which was sug- 
gested by van Zee et al. (1998b). This picture can explain the 
extremely low metallicity of I Zw 18, but is in contrast with the 
results of Aloisi et al. (2007), who concluded that I Zw 18 has 
old stars and is not a young galaxy in formation. Alternatively, 
I Zw 18 A and I Zw 18 C may be old stellar systems that 
are accreting cold gas from the inter-galactic medium and are 
now forming new stars. This may be in line with some simula- 
tions of dwarf galaxy formation (e.g. Keres et al. 2005; Dekel & 
Bimboim 2006). 

6.2. Comparison with other dwarf galaxies 

The evolution of BCDs is still an open issue. In particular, it is 
not clear what objects can be identified as their progenitors and 
descendants (Papaderos et al. 1996; van Zee et al. 2001). It is 
useful, therefore, to compare their properties with those of other 
types of dwarf galaxies. 

In Fig. 10 we compare I Zw 1 8 A with a typical dwarf irregu- 
lar (UGC 7232), taken from the sample of Swaters et al. (2009). 
The H I observations for these two objects have almost the same 
linear resolution (-200 pc), making it possible to compare HI 
surface densities and velocity gradients. The two galaxies have 
approximately the same H I size and the same rotation velocity 
at the last measured point, thus they have roughly the same dy- 
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Radius (kpc) 

Fig. 10: Comparison between I Zw 18 A (blue line) and a typical 
dwarf irregular UGC 7232 (red line), selected from the sample of 
Swaters et al. (2009). Top: HI surface density profile. Middle: 
R-band surface brightness profile. Bottom: H I rotation curve. 

namical mass. However, their structural properties are very dif- 
ferent: 

1 . the azimuthally averaged H 1 surface densities of I Zw 1 8 A 
are a factor ~2 higher than those of UGC 7232 (Fig 10, top). 
Also, the Hi distribution of I Zw 18 A is clumpy with HI 
column densities as high as -50-100 M© pc"^ (Fig. 3); 

2. the stellar component of I Zw 18 A is much more compact 
than that of UGC 7232 (Fig. 10, middle); 

3. the rotation curves have a completely diff'erent shape (Fig. 
10, bottom). UGC 7232 has a slowly rising rotation curve, 
indicating a smooth mass distribution. I Zw 18 A, instead, 
has a flat rotation curve with an inner steep rise, indicating 
the presence of a strong concentration of mass. 

In Section 4.2 we showed that the central mass concentration 
cannot be explained by the newly formed stars and/or by the con- 
centration of HI, but it can be identified with the old stars and/or 
dark matter. As to the molecules, their amount is highly uncer- 
tain and it is not clear whether they are dynamically important or 
not. The now uncovered concentration of mass is unique amoung 
dwarf irregular galaxies and must be tightly linked with the star- 



burst. Furthermore, this result sheds new light on the question of 
the evolution of BCDs and of their descendants. It is clear that, 
unless a significant re-distribution of mass takes place, a steeply 
rising rotation curve, as found for I Zw 18, would be the distinc- 
tive signature that makes the descendants recognizable. 

Regarding other BCDs, it is known that their underlying stel- 
lar component, which is formed by old stars, is generally more 
compact than common dEs and dirrs (e.g. Papaderos et al. 1996; 
Gil de Paz & Madore 2005). In particular, the old stellar com- 
ponent of the majority of BCDs has a typical central surface 
brightness /io ~ 21 mag arcsec"^ in the B-band (e.g. Gil de Paz 
& Madore 2005), similar to high surface brightness (HSB) disk 
galaxies (e.g. van der Kruit & Freeman 201 1). If the distribution 
of mass is strongly coupled to the distribution of light (Sancisi 
2004), we expect that BCDs show a dynamical behaviour sim- 
ilar to that of HSB spiral galaxies, i.e. steeply rising rotation 
curves that can be described under the maximum disk hypothe- 
sis. This seems to be the case for I Zw 18 A. In all these respects, 
I Zw 18 A resembles a "miniature" HSB disk galaxy. There are 
already indications that BCDs may have "steeper rotation curves 
than similar luminosity, low surface brightness dwarf galaxies" 
(van Zee et al. 2001), but a detailed dynamical study is needed to 
trace reliable rotation curves and determine the relative contribu- 
tions of gas, stars and dark matter to the gravitational potential. 

7. Conclusions 

We analysed H I observations of the blue compact dwarf galaxy 
I Zw 18. Our main results can be summarized as follows: 

- the H I associated with the starburst region (I Zw 1 8 A) is in 
a compact rotating disk. The H I column densities are very 
high, up to -50-100 Mo pc'^ (~ 0.6-1.2x10^2 atoms cm'^); 

- the disk has a flat rotation curve with an inner steep rise. 
This indicates the presence of a strong concentration of mass, 
that may be luminous or dark. Baryons may dominate the 
gravitational potential in the inner regions; 

- the disk has a radial inflow/outflow motion of ~15 km s"^ ; 

- the stellar concentration to the North- West (I Zw 18 C) is 
surrounded by extended H I emission, smoothly connected 
withlZw 18 A; 

- a HI tail extends to the South of I Zw 18 A out to -13.5 
kpc. It shows a coherent kinematical structure and seems to 
be connected with the H I emission to the North- West. 

I Zw 18 A appears structurally difl'erent from a typical dirr 
in terms of HI distribution, stellar distribution and dynamics. 
In particular, it has a strong central concentration of mass. 
It may be considered as "miniature" HSB disk galaxy. The 
H I concentration and the dynamical properties must be tightly 
linked with the starburst. They are also crucial to address the 
question of the progenitors/descendants of BCDs. 

Regarding the mechanism that triggered the starburst, an in- 
teraction/merger between gas-rich dwarf galaxies seems to be 
the most likely hypothesis. 
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